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1. Introduction 


Let A be an N x N matrix that satisfies the power boundedness 
condition 

p(A) = sup IIA’^11 < (1) 

n > 0 

where II *11 = ® power series expansion it Is readily verified that 

A then also satisfies the resolvent condition 

r(A) = sup ( |z |-l)ll (zI-A) ^11 < ”, (2) 

|z|>l 

and moreover r(A) < p(A). One of the assertions of the Kreiss Matrix Theorem 
[3,4,7] is that the converse is also valid: if r(A) < ”, then p(A) < “ 

also, and p(A) can be bounded In terms of N and r(A) but otherwise 

independently of A. This result is useful in proofs of stability theorems 
for finite difference approximations to partial differential equations. 

In this note we resolve an old question that has been contributed to most 
recently by Tadmor [8]: given N and r(A), how large can p(A) be? 

According to Tadmor, Kreiss's original proof in [4] unwinds to give a far from 
sharp bound 

p(A) < [rCA)]*" (VA) 

which subsequent improvements by Morton, Strang, and Miller lowered to 

2 

p(A) < 6^(N+4)^”r(A), /r(A), e^'^ r(A) (VA). 
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A few years ago Strang (private communication) observed that a paper of Laptev 
[5] implicitly derives a much more reasonable estimate [3] 

p(A) < ^^r(A) (VA). 

Finally Tadraor"s proof, which makes use of an elegant Cauchy Integral argument 
adapted ..from Laptev, yields a bound that is linear in N, 

p(A) < ^l^r(A) (VA) (3) 

Tadmor conjectures that a linear dependence as in (3) is the best possible. 
However, up to now the strongest growth of p(A) with r(A) attained by an 
example has been logarithmic, i.e., p(A) » r(A)logN [6]. 

First we will show that Tadmor's conjecture is correct, by exhibiting a 
family of matrices which p(^) ~ N+“>. By refining 

the Cauchy integral argument, we will then show that for arbitrary matrices 
(3) can be sharpened to p(A) < 2eNr(A). (Our proof is essentially Tadmor's* 

but gains the factor 16 /tt over his by dealing with complex functions 

directly rather than taking real and imaginary parts.) Together these results 
establish that eN is the optimal constant of proportionality relating 
p(A) to r(A) except for a possible factor of 2. The final section will 
prove analogous results for the continuous problem involving matrix 
exponentials e^*". 



3 


2. tiMBple with p(Ajj) ~ ®Nr(Aj^) . 

Consider the N x N Jordan matrix 




0 Y 

0 Y 


0 Y 

0 


with N > 3, Y > N (these constraints could be relaxed considerably). For 
this matrix one has IIA^II = y” for n < N-1 and IIA^II = 0 otherwise, so 


A is power bounded with 


P(A) = Y 


N-1 


(4) 


On the other hand the resolvent matrix is 


(zI-A) 


-1 


1 y/z (y/z)‘ 


1 y/z (y/z)' 


(y/z) 


N-1' 


From the fact that II Bll < upper-triangular Toeplitz matrix B, 

we obtain with a little calculation the estimates 


2^/lzl if Iz| > y/2. 


(zI-A) II < 


By (2), one therefore has 


y'^ ^(1-|z|/y) ^|z| ^ if |z| < y/2. 


r(A) < max 


sup (p-I)y^ ^(1-p/y) sup (p-l)2^/p 

1<P<Y/2 P>y/2 


“•1 “"2 

This maximum is attained at a point p - 1+N ^+0 (n ), where the estimate 


becomes 
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r(A) < 




(5) 


since Y > N. Comparing (4) and (5) shows that for this example one has 


P(\) ^ (eN-const)r(Ajj) , 

as required. 


( 6 ) 


3. Proof of p(A) < 2eNr(A) for all A. 

Theorem 1. Let A be an N x N matrix with r(A) < <». Then 


p(A) < 2eNr(A). (7) 

Remark . The factor of 2 is probably unnecessary; see the remark after 
the Lemma in the Appendix. 

Proof. Suppose r(A) < ». The matrix A” can be written in terms of 
the resolvent by means of a Cauchy integral (see [2], pp. SSS-S??) 

a" “ ^/z"(zI-A)"^dz, (8) 

where the contour of integration is any curve enclosing the eigenvalues of A, 
which must all He in lz| < 1 since r(A) < «>. Let u and v be arbitrary 
unit N-vectors, i.e., Hull = llvll = 1. Then 

* n 1 f n . . 

^ ^ q(z)dz 
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where q(z) == v*(zI-A)”^u. Integrating by parts gives 

* . n -1 f n+1 , , . j 

* " ■ jiKsmi’’ ’ 

Let the contour of integration be taken as | z | “ 1 + On this path one 

has Iz’^^^l < e, and there follows the bound 

Iv A”u| < J 1 k'(z)IIdz|. 

2.(n+l) 1 

n+1 

Now as verified on p. 155 in [8], q is a rational function of degree N. By 
the lemma in the Appendix, the integral above is accordingly bounded by 4itN 
times the supremum of lq(z)| on |z| = 1 + l/(n+l), and by (2) this supremum 
is at most (n+l)r(A). Hence we obtain 

|v A*^ul < 2eNr(A). 

Since IIA'^II is the supremum of |v*A*Si| over all unit vectors u and v, 
this proves the theorem. | 

4. Aaalogous Results for e^^ 

For problems that are continuous in time rather than discrete, stability 

At* 

depends on the boundedness of a family of matrix exponentials e (t>0) 

rather than of powers A”. Correspondingly, the resolvent of A is of 
Interest for z in the right half plane rather than outside the unit 
circle. Following (1) and (2), define 
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and 


P(A) = suplle'^*'ll 
t>0 


R(A) = sup Rezll(zl-A) ^11. 
Rez>0 


( 10 ) 


As before, one has R(A) < P(A) readily, this time as a consequence of the 

oo 

-1 r -Zt At 

Laplace transform formula (zI-A) “ J e e dt. The continuous form of the 

0 

Kreiss matrix theorem asserts that conversely, P(A) can be bounded in terms 
of R(A) and N, independently of A. We make this sharp by essentially the 
same argument used before: 


Theorem 2 . Let A be an N x N matrix with R(A) < «• Then 

P(A) < 2eNR(A). (11) 


Remark » Again the factor of 2 is probably unnecessary. The 

Laptev/Tadmor estimate has a constant 32e/Tt, as in (3). 


Proof. In analogy to (8), one has now 


e^*^ = :^/e^’^(zI-A)"^dz, 


where the contour of integration can be taken as any line Rez “s' p>0. 
Integration by parts gives 


* At 
V e u 


1 . 


2irit 


-/e q" (z)dz 


with q(z) = v*(zI-A)~^u. Taking the contour U = now leads to the desired 
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bound (11), again by the use of the lemma in the Appendix. | 

Constructing an example to prove that (11) is sharp, on the other hand, 
is trickier than it was in the power-boundedness case. The following example 
achieves growth proportional to /N, not N. Omitting details, define 



Then one has 


N-1 


yt 


(yt) 
(N-l)! 


At -t 
e = e 


yt 


L ^ J 

For large y , this matrix achieves maximum norm near t = N, where it is 
dominated by the upper-right entry, with magnitude approximately 




e 


-NjjN-l^N-1 

~(N-1)! 


,N-1 


✓ 2tiN 


( 12 ) 


For the second estimate we have used Stirling's formula. On the other hand 
the resolvent matrix is 


(zI-A) 


-1 


z + 1 


z + 1 
1 


z + 1 


''Z+1-' 


1 
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For large Y, Rez times the norm of this Is maximized near z = 1/N, where 
again the upper-right entry dominates and one has 




N-1 


eN 


(13) 


Comparing (12) and (13) shows that in this example one has 


Appendix — Lemma on arc length of a rational function on a circle. 

Let S be any circle or line in the complex plane, and define the 

and norms over S by Il f II = / If(z)||dzl, II fll^ = suplf (z) I . The 

S S 

following lemma provided the key argument in proving Theorems 1 and 2. For 
the case of a polynomial the result is a corollary of Bernstein's inequality 
[1], llq^ll^ < Nllqll^^ for S = {z:Iz|=l}, but the extension to rational 

functions appears to be new. Since Hq'll^ represents the arc length of the 
image of S under q, the lemma has a simple geometric meaning. The 
example q(z) = b(z-s), where b(z) is any finite Blaschke product of degree 
N (such as z~^) and s is the center of S, shows that it is sharp except 
for a factor of 2. 


Lemma . Let q be a rational function of degree N with no poles on 
S . Then 


llq'll^ < 4iTNIIqll^. 


Remark. We believe that the bound is valid with a factor 2ir instead of 


4tt, but have been unable to prove this. 
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Proof ♦ Since the composition of q with a Mobius transformation is 
again a rational function of type N, we can assume without loss of generality 
that S is the unit circle. Define g(z) to be the angle of the tangent 
to q(S) at q(z), l.e. 


g(z) = arg[zq^(z)] . 

Let TV[g] be the total variation of g over S. The lemma is a consequence 
of the following two facts: 

(a) llq-ll^ < TV[g]llqll^, 

(b) TV[g] < 4ttN. 

The proof of (a) is a matter of integration by parts: 
llq^ll ^ ^ Iq'(z) I Idz| = jzq'(z)e 

= -ifq'(z)e = ^q(z)g' (z)e ^®^^^dz 

< llqll Jlg"(z)| Idzl = llqll^TV[g]. 


To prove (b), note that q' is of rational type (2N-1,2N), so zq^(z) 
is of rational type (2N,2N) and can be written as a product 


zq'(z) 


2N \ 


z + b. 


k=l 


c, z + d, 
k k 


g(z) 


2N 


I arg( 
k=l 




This Implies 
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and therefore 


TV[g] < 


2N 

I TV[arg( 
k=l 




.)] < 4ttN. 
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